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I. INTRODUCTION 



Supernoave type la data [T] as well as other observational probes [2] show that the Universe is undergoing an 
accelerated phase of expansion at present time, a feature that does not emerge from standard cold dark matter 
model [3 . Attempts to explain such unexpected behaviour go towards modifications of either the geometric part 
of the Einstein field equations, implying modified theories of gravity, or the matter sector, thus involving new and 
sometimes weird forms of energy Such a dark sector seems unavoidable in order to fit present cosmological data. 
On the other hand, from the theoretical side, the strong mathematical resemblance between the dynamics of spacetime 
horizons and thermodynamics is strongly attested [S] |B] so that gravitational fields equations can be given a physical 
interpretation which is thermodynamical in origin. In particular, the Friedmann equations follow from applying 
the first law to the apparent horizon of an isotropic and homogeneous universe, not only in Einstein gravity, but 
also in more general Lovelock gravity [7]. Likewise, it seems that a gravitational theory built on the principle 
of equivalence must be thought of as a macroscopic limit of some underlying microscopic theory, the microscopic 
structure of spacetime manifesting itself only at Planck scale or near singularities. Also the horizons link some aspects 
of microscopic physics with the bulk dynamics [5]. It is well known since long that one can define entropy and 
temperature for a spacetime horizon [SHU]; in fact, many attempts have been done to better understand this link. An 
instructive example is the case of spherically symmetric horizons in four dimensions, for which Einstein's equations 
can be interpreted as a thermodynamic relation arising out of virtual displacement of the horizon [12j . Moreover, the 
same interpretation holds for the case of the Lanczos-Lovelock gravitational theory in D dimensions [13] and explicit 
demonstration has been given for Friedmann models |14l. 

In the present paper, following this deep relation between thermodynamic and gravity, and in particular between 
entropy and horizons, we argue that some form of dark energy is demanded on thermodynamic grounds. 
In order for an isolated system to evolve to thermodynamical equilibrium, the entropy function of the system must 
show two properties: it must never decrease, i.e. its first derivative with respect to the relevant variable must be 
non- negative, and convex, i.e., its second derivative must be negative. 

This constitutes the hard core of the second law of thermodynamics and it is naturally realised in systems dominated 
by electromagnetic forces; however it might not be true when gravity plays a role. In fact, the entropy of the system 
must still increase but it may be grow unbounded: this occurs, in the Newtonian framework, for the Antonov's sphere, 
the final stage of N gravitating point masses enclosed in a perfectly reflecting, rigid, sphere whose radius exceeds some 
critical value [151 [TB]. Nevertheless, when we replace Newtonian gravity by general relativity, a black hole is expected 
to be formed at the center of the sphere that, though it tends to evaporate, it will likely arrive to an equilibrium state 
characterised by a state of maximum, finite, entropy. 

In any case, in a Friedmann- Robertson- Walker (FRW) cosmology, the Universe seems to behave as an ordinary system 
whose entropy increases towards a maximum value. The latter follows from the observational data on the evolution of 
the Hubble factor of the FRW metric [T7] and from the evolution of the entropy of the apparent horizon, that seems 
to be the appropriate thermodynamic boundary [IB] . 

The present paper is a second step of the analysis outlined above. Infact, in a previous paper |19j we showed that 
an Einstein Universe, as a thermodynamical system, cannot tend to equilibrium in the last stage of expansion unless 
it accelerates. We have found that this holds true for some modified models that are dynamically equivalent at the 
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background level, nevertheless this does not mean that every accelerating universe is thermodynamically motivated 

In this work we study the thcrmodynamical behaviour of a braneworld model with two correction terms: a four- 
dimensional curvature on the brane and a Gauss-Bonnet (GB) term in the bulk |20j . The induced gravity (IG) 
correction arises because the localized matter fields on the brane, which couple to bulk gravitons, can generate via 
quantum loops a localized four-dimensional world- volume kinetic term for gravitons [21j . On the other hand, a Gauss- 
Bonnet term naturally appears in an effective action approach to string theory, corresponding to the leading order 
quantum corrections to gravity [22 . As a result, we have the most general action with second-order field equations 
in five dimensions [53]. 

Section |ll] introduces the braneworld cosmology with induced gravity and curvature effects. Subsections |II B| and 
|II C| focus on the entropy of the horizon and matter components, respectively. The energy componentsof the Universe 
are assumed to enter the field equations in the form of perfect fluids, the standard equation of state being true for 
each of them: pi = WiPi. Then, in section III the matter component is assumed as cold matter and a Chaplygin 
gas. The choice of a Chaplygin gas is based on the recent observational result that the equation of state parameter 
of dark energy can be less than —1 and even display a transient behaviour [5]. This can be achieved either by 
means of phantom fields, that on the other hand suffer from instabilities |24j . or by other approaches that mimic this 
phantom-like behaviour. In the model under analysis, in which UV modifications are included by considering the 
stringy effect via the GB term in the bulk, and IR modifications are due to the IG effect, a Chaplygin gas fluid on the 
brane provides a smooth crossing of the cosmological constant line j25| . In fact, this component is characterised by 
a cross-over length scale below which the gas behaves as pressureless fluid and above which it mimics a cosmological 
constant. 

Our conclusions are given in section |IV[ we find that even this modified theory of gravity needs a component with 
typical dark energy behaviour in order to satisfy the generalised second law (GSL) and approach thcrmodynamical 
equilibrium in the long run. 



II. GAUSS-BONNET AND INDUCED GRAVITY CORRECTIONS 

The total action of the braneworld model under consideration reads [ID] 

/ = ^ y rf'a; (^"^R - 2A5 -f aCoB) + ^ / {R - 2A4) + j d^x C^n, (1) 

where A5 < is the cosmological constant on the bulk and 

Cgb = ^'^R' - 4 (^'i?^^ ^"^Rab +('^ i?^^^^ ^"^Rabcd (2) 

is the GB correction term, whose coupling constant a = l/Sg^ is related to the string energy scale, gs [IHj- The 
gravitational coupling constants K4 = 87rG'4 and k| = SttG^ on the brane and in the bulk, respectively introduce a 
length scale, the induced gravity cross-over scale, r — n'l/2n\ and help defining the brane tension, A = A4/K4. Last 
term represents matter action. 



A. Cosmological equations 



The metric of homogeneous and isotropic FRW Universe on the brane, with spatial curvature index fc, is 

ds^ = h^^dx'^dx" + [d9^ + smO^dct)'^] , (3) 



where f = a{t)r, the two-dimensional metric is h^^ — diag (— l,a^/(f — fcr^)) with x'^ — t and x^ — r. This allows 
the explicit evaluation of the radius of the apparent horizon (a marginally trapped surface with vanishing expansion) 
determined by the relation h^^^d^rdi,f = [SB] that gives 
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Friedmann's equation on the brane is 
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in which, assuming there is no black hole in the bulk, — j;^ (^^^ + y^l + f^As^ and matter field are supposed 

to be perfect fluids with energy density pi so that in order to describe completely the cosmological dynamics on the 
brane we can use the energy conservation law 



where Wi — Pi/pi = const. 



Pi + 3Hp,{l + w,) ^0, 



B. Entropy of the apparent horizon 
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The entropy on the apparent horizon is given by 
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Its derivatives with respect to the scale factor a, that will be denoted by a prime, can be calculated and simplified by 
using eqs. 
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where wt = Pt/pt = J2tPi/J2iPi- 

It can be easily checked that the entropy grows provided the fluids component of the Universe satisfy wt > — 1- 
In order to evaluate the second derivative we make use of the late time behaviour of this brane cosmology with 
curvature correction that, as found in '205 1 reduces to conventional cosmology with positive effective gravitational and 
cosmological constants. For a — > cx) cq. ([s]) reduces to 
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/3, 7 and v being functions of the parameters of the model, namely a, A4, A5 and r, - see |20j . By using this solution 
the second derivative of the horizon entropy w.r.t. the scale factor, it gives 
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where wi is the equation of state parameter of the fluid that redshifts more slowly. Thus, entropy tends to equilibrium 
if Wt > — 1, which is required by the GSL. 



C. Entropy of matter fields 



We first consider a perfect fluid with Wi ^ 0. In this case, using Gibbs' law [30] 



TfdSf = d{p,V)+p,dV, 



(13) 



^/o. The evolution of the fluid temperature is governed by dlnTf /dlna — —3wf, see e.g. |31j . whence 
3wi ^]-^jg jxieans that fluids with Wf > dilute as the universe expands and the temperature decreases 



with V 
Tf = Tfoa 

while dark energy fields behave the opposite and the temperature will grow during the expansion; for details see 
We then compute the two first derivatives of the entropy of the fluid enclosed by the apparent horizon 
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A simple evaluation of the first derivative for a — )■ cxd shows that S'f 



^ when 1 + > 0: it tends to zero from 



below. To discern which component prevails in the long run we evaluate the ratio between the derivatives: 
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provided — 1 < ui; < 0. The same holds true for the ratio Sj/S'J^. 

Last section is devoted to dust fluid, for which we cannot use the evolution of temperature in terms of the scale 
factor as above because it gives the unphysical result of a vanishing constant temperature. Then, we proceed as in 
[TO] in which every dust particle contributes a given bit to the fluid entropy so that Sm = N with N = 47rf^n/3 
being the number of particles within the apparent horizon and n = noa^'^ the number density of dust particles. 
By computing first and second derivatives of dust entropy inside the apparent horizon we find 
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Now, by computing the ratio between the components as before, we obtain 
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It means that if the only fiuid component is dust, we need to compute explicitly the sum of the horizon and matter 
contribution to the entropy derivatives. By restoring proper dimensions we get 
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In order to fulfill the GSL the right hand side of eq.(19l must be always positive and to approach the equilibrium 



eq. ( 20 1 must be negative in the last stage of expansion. These two requirements translate in the following constraint 
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Now, since pmo = n^rnc , where m is the mass of the particles of the cosmological fluid, the apparent horizon radius 
must be larger than the associated reduced Compton wavelength of the particles, Ac = fi/mc. 

If, on the contrary, we had more components, > would imply that the dust term would dominate in the long 
run and its derivatives fulfill neither the GSL nor the equilibrium prescription. By requiring — 1 < < the horizon 
entropy will prevail in the long run and the Universe tend to equilibrium, as the analysis of the previous section shows. 

To sum up, a dark energy component (i.e., a barotropic fiuid with —1 < w < 0) in the late time evolution is still 
needed, in such a modified theory of gravity, if the Universe is to tend to equilibrium eventually. 



III. A DARK ENERGY COMPONENT: CHAPLYGIN GAS 



As concluded is the previous section, it seems that from a thermodynamical point of view a dark energy component 
is needed. In what follows we restrict ourselves to a Minkowski bulk and a spatially fiat FRW brane so that the radius 
of the apparent horizon now reduces to ta — l/H. When A5 = 0, $0 could be zero or $0 = — 1/2q;; we use the first 
choice since the latter case has no IG limit: it corresponds to an AdS bulk even though A5 = [33 . Cosmological 
evolution still follows eqs. (l5])-(|6]) where the total energy density is now 
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Observations 




As 






SNIa (gold sample) 


0.43 0.11 


0.50 


1.45 0.99 0.99 


0.923 


SNIa(gold sample)+CMB+SDSS 


0.51 0.41 


0.099 


1.79 0.99 -0.98 


0.992 



Table I: Best fit values for the parameters of the GBIG Chaplygin model by fitting with SNIa gold sample and 
SNIa+CMB+SDSS. Cosmological parameters are defined as VLr = l/4rff^, Q.^ = 8q//^/3, ^a = k1A/3//^, and ^^h = 
PchoK-l/^Ho. fi™, = kIA/3Ho has been fixed to Qm = 0.0456, according to WMAP and BBN resuhs |25] . 

and the CDM and Chaplygin gas components obey: 

p,n + 3Hp„, - 0, (23) 
Pah + SH{pah+Pch) ^0, (24) 

respectively. Then CDM redshifts as p„i ~ PmViar^ while, being the Chaplygin pressure pch = WchPch — ~-^P^h i 
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where As = ApJ^'^"^ , 



and alpha a constant. 



A. Entropy of the apparent horizon 



The first and second derivatives of the apparent horizon entropy Sa for such a model read 
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Let us evaluate the sign of these derivatives: in order to fulfill the GSL the sum of the energy density and pressures 
of CDM and Chaplygin must be positive. To estimate it we use the best fit results obtained in [25 whose model 
was tested using data from SNIa (gold sample), CMB (shift parameter) and SDSS (BAO peak). Qualitative different 
evolution and best-fit results are obtained for SNIa data and the joint analysis of the total data set, as can be seen 
in Fig[l]and table [l] the two panels plot the evolution of the energy densities and pressure multiplied by 3Hq / {SttG): 
Pm, Pch{^ + Wch), and the sum of the two. On the left panel best fit values for the solely SNIa fit has been used 
and a slower redshift of the CDM can be appreciated. On the contrary, when SNIa-|-CMB+SDSS fit is taken into 
account, as displayed in the right panel, the Chaplygin gas redshifts more slowly. Then, in what follows we will test 
our thermodynamic hypothesis using both fits. 

Moreover, Figjl] shows that the term in brackets is always positive which ensures the validity of the GSL. 
The second derivative is negative, the minus sign given by the terms in squared brackets that, for SNIa fit tend to 
—4, while in the other case tend to —4 — 3a ~ —1.06: the term with the derivative of the horizon radius tends to zero, 
since the model tends to De Sitter, and the term involving energy densities tends to zero for the SNIa fitting and to 
-35 for the SNIa-^CMB-^SDSS one. 



B. Entropy of matter fields 



In this section we evaluate the role of the CDM entropy. We note that eqs.(17)-(18l still apply. 
Computing, as before, the ratio between CDM and horizon entropy derivatives one can see that for the best-fit values 
obtained just using the SNIa data 

^ ^ ^ ^ const (28) 



SO that eqs.(19)-(20l as well as the constraint below still hold. 

On the other hand, when computing the ratio in the second fit case, due to the fact that Chaplygin gas component 
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Figure 1: Qualitative evolution of the energy densities and pressure pm (dot-dashed line), pch + Pch (dashed) and the sum of 
the two (solid) for different values of the parameters of the GBIG Chaplygin gas model. The left panel shows the evolution 
for the best-fit values by fitting the model with data from SNIA (gold sample). In the right panel all the observational data 
(SNIa, CMB and SDSS) are used in the fit. 



redshifts slower than CDM, both first and second derivatives ratios tend to zero. In this case the horizon contribution 
prevails in the long run ensuring that both the GLS is fulfilled and the equilibrium will be eventually approached. 
As before, we can argue that only in the case that a dark energy component eventually dominates, the second law of 
thermodynamics is fulfilled. 

Finally, we evaluate the Chaplygin gas contribution to the entropy. It yields 
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The temperature evolution is found by solving = —3wchTch/a. It turns to be 



Tch = Tcho (l - A, + ^a^(i+": 



1/(1+5) 



(31) 



For both best-fit we find that 



qi c/r 

and 



(32) 



when a — i' oo. Thus, this contribution does not change previous results: a dark energy contribution must be present 
in the model and must dominate in the long run if the Universe is to eventually tend to thermodynamical equilibrium. 



IV. CONCLUSIONS 



We have shown that a Gauss-Bonnet model with induced gravity corrections tends to thermodynamical equilibrium 
in the long run provided that at least one of the cosmological fluids corresponds to dark energy. The apparent horizon 
entropy is well-behaved, i.e. it is a growing and convex function of the scale factor, while the entropy of the fluids, 
either cold matter or radiation, is growing but diverging. Nevertheless, if a dark energy fluid {w < 0) prevails in the 
long run, the overall entropy tends to a finite maximum value. 

An explicit example is given by a Chaplygin gas component that behaves as a dynamical dark energy, smoothly 
crossing the phantom divide line and approaching the cosmological constant behaviour. The model under analysis is 
dominated by two energy components, CDM and Chaplygin gas. Best fit values for the parameters of the model were 
obtained in [3S] and very different results are found when testing the model against SNIa data or SNIa-|-CMB-|-SDSS 
set. This implicitly translates into a different weight of the Chaplygin gas component in the long run thus altering the 
results on the entropy functions. When we use the best-fit values of the SNIa data alone, cold matter entropy evolves 
the same way as the horizon entropy, then in order to fulfill thermodynamical requirements the horizon radius must 
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be larger than the Compton length of the cold matter particles. On the other hand, when we use SNIa+CMB+SDSS 
best fit results, since the Chaplygin gas redshifts more slowly, the overall entropy in the long run reduces to the 
horizon entropy that evolves toward equilibrium. 
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